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I.  INTRODUCTION 

In  this  sLudy,  we  sre  concerned  vit.i  Che  determinetiou  end  grephlcsl 

presentation  of  a  self-coneietcnt  see  of  blast  wave  properties  for  a  TNT  exo’-'Sive. 

These  date  can  then  be  utilised  to  predict  the  characteristics  of  the  resulting 

blast  field  .nvironnent.  It  is  to  be  noted  that  a  similar  study  by  Baker  and 
1* 

Schuman  can  b«  ueed  to  predict  the  blest  environment  generated  by  a  Pcntollte 
exploalve  source. 

The  following  sections  of  this  report  will  review  the  available  sources 
of  experimental  data  on  bleata  from  both  TNT  and  Pentolite  in  order  to  provide 
the  reader  with  a  comprehensive  listing  of  such  sources  of  date.  The  TNT  data 
from  a  selected  nua&er  of  these  references  will  be  coshered  to  the  various 
theoretical  prediction  techniques  that  ere  available,  in  order  to  define  the 
fector(s)  necessary  to  effect  agreement  between  the  experimental  date  end  the 
theoretically  predicted  magnitudes  of  these  date.  Finally,  e  series  of  graphs 
will  be  presented  that  will  provide  e  self-consistent  set  of  air  blest  parameters 
for  explosive  datoneM".,*  of  TNT  In  air.  Uhera  appropriate,  methods  of  account¬ 
ing  for  ground  reflection  effects  sod  conversion  from  TNT  to  Pentolite  explosive 
will  be  discussed. 

The  previous  study  conducted  by  Baker  and  Schuman^  resulted  in  e  s'.itler 
set  of  elr  blast  parameters  for  Pentolite  detonellons.  However,  th»s  study  wes 
based  on  correlating  the  shock-front  parameters,  es  compiled  by  Goodman*  from 
experimental  date  obtained  from  Pentolite  detonations,  with  plots  o.  date  behind 
the  shock  front  obteined  from  the  theoretical  predictions  of  Brode*  tor  TNT 


*  Superscript  numbers  denote  references  in  Section  VII. 
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detonations.  Adjustments  were  made  hy  the  author*  ot  reterenca  1  to  account  tor 
the  d  'ference  In  type  of  explosive*,  and  to  match  the  theoretical  analyse*  with 
the  experimental  date  along  the  shock  front.  Hteever,  the  theoretical  TNT 
analytes  could  not  be  adjusted  to  exact  coincidence  with  the  PcnLolite  expert' 
mental  results;  furthermore,  the  analytes  of  Brode^  exhibited  Internal  inccr.ais- 
tenclca  that  could  not  be  overcome. 

Tn  this  study,  we  will  attempt  to  eliminate  some  of  the  Inconsistencies 
of  reference  1  by  comparing  TNT  experimental  date  with  e  theoretical  prediction 
technique  baaed  on  e  YNT  detonation.  It  It  to  be  noted  that  the  theoretical 
analysis  selected  In  thle  study  doss  not  suffer  from  Internal  Inconsistencies 
ee  doee  the  data  predicted  by  Brod*\  Moreover,  thle  theoretical  approach 
appears  to  be  the  bast  available  at  the  preeent  time  In  that  It  provides  batter 
agreement  with  actual  experimental  date  than  doee  the  theory  advanced  by  Brode. 

The  self- cons  latent  eat  of  bleat  parameter*  presented  In  this  report 
have  bean  generated  for  MSA  under  MSA  Contract  No.  bA*49>  146-VXZ-29J.  In 
addition,  the  theoretical  computation*  were  perforated  by  NOL  using  their 
‘Vundy  75"  Cod*.  AAI  wishes  to  ecknot  ledge  the  ki-.~  coop- ration  of  personnel 
of  both  of  these  agencies  during  preparation  of  this  report. 


In  Pact  A  of  this  section  we  will  surmarlse  the  content*  of 
various  references  chat  provide  experimental  data  defining  the  blast 
wave  characteristics  that  result  from  detonations  of  both  Pentoll'te  an : 

TNT  explosives.  The  data  obtained  from  a  selected  number  of  the  TNT 
references  will  be  used  In  latte”  •  tl  or.s  of  this  report  to  graphically 
formulate  the  self-consistent  set  ot  blast  wave  properties  for  this 
explosl«*. . 

It  Is  to  be  noted  that  the  data  sources  discussed  In  this  section 
do  not  represent  an  exhaustive  listing  of  all  available  references; 
however,  these  sources  have  been  judiciously  chosen  to  provide  a  representative 
cross-section  of  the  presently  available,  •r-erlmentally  derived.  Information. 
These  data  sources  are  referenced  In  Section  VII  of  this  report. 

Following  the  discussions  of  there  data  references,  a  compre¬ 
hensive  chart  that  deflnas  the  type  of  explosive,  ambient  condltlais,  the 
scaled  distance  range  covered  by  the  collected  Hats,  and  the  type  of 
blast  wave  properties  measured  in  each  rvfvrencc  discussed  In  this  section, 
is  presented. 

In  Part  I  of  this  section,  an  evaluation  of  these  data  source 
references  is  given  In  terms  of  the  applicability  of  the  Information 
contained  in  these  i*f«iences  to  the  forauletlon  of  a  self-consistent 
set  of  blast  wavs  properties,  The  conclusions  drawn  In  this  discussion 
will  dictate  the  selection  of  dara  r>>  he  used  in  the  latter  sections  of 
this  report. 
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A.  Data  Source  Discussion 

The  graphical  airblast  data  presented  by  Baker  and  Schuman' 
was  developed  primarily  for  use  in  predicting  free  flow  blast  conditions 
on  moving  airfoils.  The  shock  front  parameters  presented  in  reference  1 
have  b»en  extracted  iron;  the  compiled  experimental  Pentolito  data  given 
by  Goodman^,  whereas  the  data  for  parameters  behind  the  shock  front  were 
obtained  from  the  theoretical  predictions  of  Brode^.  Unfortunately. 

Brode's  data  could  not  be  adjuated  to  provide  exact  agreement  with  the 
experimental  data.  P.eference  1  recommends  the  following  simple  conversion 
factors  to  account  for  variation  in  typr  of  explosive  and  tor  effect  of 
ground  reflection: 

1.1  lb  TNT  as  1  lb  SO/ 50  Pentollte* 

and  Ftee  Atr  Weight  of  Pentol Ite  «  1  .*  ••entolite  weight 

on  the  ground 

The  compiled  experimental  data  reported  by  Condition  .»  based  on  th>-  av»iiabl 
air  shock  data  tor  the  1945*1960  period.  Ax  noted  In  Table  1,  these  data 
cover  a  range  uf  scaled  distances  from  0.137J  (the  charge  surface)  to 

i/3  l 

b6  ft' lb  .  (The  range  of  bla;>t  parameters  presented  by  Baker  covers 
a  scaled  distance  range  from  ?  to  60  ft/lb^\  since  this  was  felt  to 
be  adequate  for  the  airfoil  vulnerability  problem  of  primary  Interest  In 
this  reference).  The  Pentollte  explosive  density  used  by  Goodman  was 

*  50/50  Pentoiire  ■  50?.  TNT  ano  50?  PfclN  by  w.-ight 
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1.65  gm/cc.  Tlie  slde-on  Impulse  end  duration  measurements  are  less 
satisfactory  chan  'he  peek  pressure  and  reflected  impulse  measurements. 
Goodman  also  develops  analytic  expressions  for  the  side-on  and  reflected 
pressure,  and  the  slde-on  and  reflected  Impulse  variations  as  functions 
of  the  scaled  distance  by  a  least  squares  polynomlnal  fit  to  the  compi l> d 
data. 

The  data  reported  hy  Weibull  are  primarily  measurements 
of  the  arrival  timet,  shock  front  velocities,  and  positive  pressure  phase 
duration  rimes  obtained  from  TNT  explosions,  over  a  range  of  ambient 
pressures  from  approximately  4  atm  ,.o  a  vacuum  ui  a  few  mo  Hg. 

Included  <n  this  reference  is  a  comparison  of  the  shock  arrival  times 
obtained  from  different  mixtures  of  TNT  and  A1  (101,  20' ,  and  30X  Al). 

1007.  TNT,  and  PETN  e  . plosives.  Weibull  also  provides  a  set  of  formulae 
that  may  be  used  to  reproduce  the  observed  arrival  and  duration  times 
with  great  precision.  A  comparison  of  the  TNT  arrival  tliae  data  observed 
by  Weibull  with  data  compiled  at  Aberdeen  Proving  Grounds,  Ballistic 
Research  Laboratories  for  Pentnlire  explosive  Is  given  In  the  Appendices 
of  this  report. 

Doll  and  Salarn^  present  data  ohtalned  from  a  carefully 
planned  series  of  tests  using  2500  lb  of  TNT  as  the  explosive  source. 
Unfortuanrely,  the  response  tiaw  (3l5  cps)  of  the  recording  galvanometers 
appears  to  seriously  limit  the  usefulness  (and  correctness)  of  these  data. 
When  plotted  on  an  overpressure-scaled  distance  diagram  ..  data  of  Doll 
and  Salmon  are  In  disagreeiaent  with  similar  data  from  other  source:;  tin- 
peak  pressures  are  too  low,  indicating  (hat  the  recording  system  missel  the 
true  peak  because  of  an  excessively  long  rcruonsc  time  chare,  uii.tn. 
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As  noted  In  Che  Introductory  paragraphs  of  the  report  by 
Fisher  and  Pittman^,  the  primary  intent  of  thla  work  was  to  demonstrate 
the  difficulty  of  detonating  small  charges  of  TNT  and  to  describe  a 
number  of  ways  of  recognising  poor  TNT  detonation  processes.  In  addition, 
this  report  Indicates  that  placing  the  booster  half  in  and  half  out  of  the 
surface  of  a  one-pound  cast  TNT  charge,  in  contrast  to  central  Initiation, 
results  in  a  detonation  that  produces  the  accepted  peak  air  blast  pressure 
and  Impulse  magnitudes. 

Although  the  major  emphasis  of  the  data  reported  by  Baker 
et  al7  was  to  demonstrate  the  effect  of  fog  and  rain  on  the  air  blast 
characteristics  from  1000  lb  Pentolite  charges,  there  were  several  recorded 
test  events  accomplished  during  clear  weather.  The  reproducibility  of 
these  data  Is  quite  good;  from  a  series  of  eight  tests,  with  three  gages 
at  each  radial  location,  the  standard  deviations  of  the  peak  overpressures 
vary  from  0.08  to  0.55  psi.  for  a  range  of  overpressure  of  1.22  to  11.03 
psl.  A  significant  result  of  this  test  series  is  that  the  ground  reflection 
factors  computed  from  a  comparison  of  these  data  with  the  free  air  data 
of  Goodman^  cover  a  ranee  of  1.88  to  2.3 A  for  increasing  scaled  distance. 
Ovbiously,  for  a  detonation  very  close  to  the  ground,  a  reflection  factor 
greater  than  2.0  is  impossible,  since  values  equal  to  2.0  Imply  perfect 
reflection.*  These  results  indicate  that  Goodman’s  data  it  probably  in 
error  for  the  larger  scaled  distance  ranger,  since  the  data  cl  tented  by 
Baker  in  this  report  *«,  as  noted  above,  almost  free  irom  experimental 
Inaccuracies . 

*  For  detonations  above  the  ground  where  the  sealed  height  is  an  appreciable 

fraction  of  the  scaled  distance,  ic-ficclice  f  fie  I»nrs  greater  7  n 

can  be  realised. 
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Reference  S  contain*  chare*  and  graph*  * how In*  th«  affect* 
of  explosion*  at  function,  of  such  variabla*  a*  eharg*  weight  and  composition, 
type  of  weapon,  distance  from  detonation,  etc.  The  data  sheet*  contained 
In  this  report  are  designed  for  >»*#  in  prediction  of  the  vsrlous  character¬ 
istic*  associated  with  detonation  of  explosive  charges.  As  such,  this 
report  Is  not  a  data  reference  report,  but  has  been  Included  In  this 
discussion  because  of  the  unique  manner  In  »d»lch  th#  prediction  charts 
ore  designed. 

The  data  of  Koffiaan  and  Mills9  presents  msasursnsnts  taken 
from  detonation  of  50/50  Pentollte  spheres  of  l/it  -  b  lb  In  weight.  The 
si  dr-on  and  reflected  overpressure  data  and  Impulse  and  duration  time 
Mssureewnts  given  la  this  report  have  been  Incorporated  In  the  complied 
d<tta  report  by  Goodman2. 

The  reflected  lapulse  neasureaents  reported  by  Olson  were 
generated  froa  detonation  of  l  lb  Pentollte  spheres  under  reduced  aabler.t 
pressures  simulating  altitude*  up  to  100.000  feet.  Analysis  of  these 
data  Indicates  that  th*  Impulse  magnitude*  scale  according  to  Sachs'  La*/*, 
as  discussed  by  Sachs11  and  Dewey  and  Sperrassa12,  although  this  Is  to  be 
viewed  with  caution,  as  noted  by  Olson.  This  Is  primarily  due  to  th* 
fact  that  one  of  th*  i...  sis*  measurements,  when  scaled  according  to 

1/3 

Sachs'  Law  fails  Inside  the  charge  surface  (2  *  ft/ib  ).  The 

neasureaents  taken  by  Dewey  and  Sperrae**12.  and  subsequent  discussion  of 
the  scaling  procedure*  verity  the  basic  assumption*  applied  by  Sachs  In 
formulating  th*  scaling  law*.  (See  Part  A  of  th*  Appendix  section  «f  this 
report  for  a  brief  discussion  of  th*  relevant  <»last  wave  parameters  formulated 
la  terne  >i  the  Sachs'  scaling  relationships). 
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The  ci  tpilstlon  uf  blest  data  presented  by  Moulton 
incorp'.'i-ai ms  both  HE  (TNT  and  Pentollte)  and  nuclear  blast  phenomena. 

The  HE  data  presented  In  this  reference  have  been  extracted  from  various 

4 

sources;  ‘;he  data  of  Weibull  is  used  to  Illustrate  overpressure  functions; 

Fisher's  and  Pittman'  results  have  been  used  to  Illustrate  the  Impulse 

functions  resulting  ftom  detonation  of  TNT  charges,  and  the  data  of 
.  14 

Grins trom  has  been  used  to  illustrate  shock  duration  time  parameters  and 

characteristics  of  the  negative  pressure  region.  This  extensive,  three 
volume  report  Is,  at  the  present  time,  undergoing  revision  and  updating. 

The  report  by  Groves^  describes  a  technique  of  determining 
the  blast  Induced  overpressures  from  distance-time  observations  of  the 
motion  of  the  shock  front.  Included  with  the  overpressure  data  of  this 
report  is  an  estimate  of  the  percent  error  of  the  various  measurements. 

Over  a  range  of  peak  overpressures  from  approximately  3  psl  to  130  psl, 
the  percent  error  associated  with  this  set  of  measurements  Is  less  than  3Z. 

The  report  by  Dewey^  describes  a  method  of  determining  the 
particle  velocities  associated  with  a  blast  event  that  employs  high  speed 
photography  of  rsmke  trails  formed  close  to  the  charge  just  be *>'«■*■ 
detonation.  The  displacement  with  time  of  the  smoke  uaile  ran  be  measured 
from  ehe  film  record  snd  the  velocity  of  the  awoke  calculated.  Analysis 
of  these  results  gives  the  particle  velocity 'time  hletory.  The  two  5  ton 
TNT  denotations  used  In  these  tests  were  also  analyzed  to  provide  reflection 
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factors  of  1.70  and  1.88.  However,  the  tame  test  outs  is  analyzed  by 
Groves^  indicates  reflection  coefficients  of  1.76  anJ  1.60  respectively. 
Unfortuantely,  these  two  references  do  not  clearly  describe  the  data 
reduction  techniques  used  to  achieve  the  quoted  reflection  coefficient 
magnitudes,  hence,  no  Judgement  can  be  mad'  r^gardln^  the  cortectness  of 
cither  report.  Groves  compares  the  measured  IN'!  data  -Hth  tbo  theoretical 
predictions  of  Brode^  in  terms  of  reflection  coefficients,  '.ilth  the  result 
thst  Brode's  overpressure  predictions  appear  low  for  large  scaled  distances. 

A  similar  conclusion  is  reached  by  Baker ^  when  comparing  the  Fentollte 

2  3 

data  of  Goodman  with  the  theoretical  predictlone  of  Brode  . 

The  100  ton  TNT  detonation  test  discussed  in  reference  IB 
Involved  the  use  of  a  hemispherical  charge  placed  on  the  ground  surface 
and  was  designed  to  provide  measurements  of  surface  and  eub*surface  blast 
phenomena  and  to  study  effecte  of  blast  on  various  targets  and  la  tunnels 
snd  underground  chambers.  The  20  ton  TNI  detonation  test  results  presented 
in  reference  19  were  designed  to  resolve  small,  but  consistent  discrepancies 
that  appeared  in  previous  tests  and  data  collection  programs.  Complete 
discussions  of  the  various  types  of  instrumentation  used  io  given  in  addition 
M  tabular  and  graphical  forms  of  data  presentation.  The  suasuirv  report 
by  Dewey***  is  not  very  useful,  since  this  reference  was  merely  oculgned 
for  presentation  at  a  technical  meeting.  However,  this  report  Indicates 
that  peak  air  velocities  determined  by  the  smoke  trail  technique  for  B  to  too  lb 
spherical  charges  of  TNT  agree  quite  well  with  the  predicted,  or  expected 
velocities.  The  duration  of  the  positive  pressure  phase  is  greater  than 
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expected  (<-'~n«red  to  Brode’s3  theoretical  enelyeee).  The  discussion  of 
a  (3  ray  &l  .itlon  technique  used  to  determine  the  time  history  of  the 
air  density  in  the  blast  field  is  perhaps  the  most  Interesting  part  of 
this  reference. 

Although  the  report  of  Shear  and  Wright23  is  primarily 

theoretical  in  approach,  it  does  present  experimental-theoretical  correlation 

data  (for  peak  overpressure)  based  on  the  experimental  data  of  Coodaan2, 

Fisher  and  flttmsn6,  and  reference  19.  The  two  reports  by  Rudlln22,23 

contain  discussions  of  the  basic  phyalcs  of  the  explosive  prodess  during 

the  early  stages  of  the  detonation  wave  and  subsequent  airshock  phase. 

The  photographic  data  obtained  for  the  very  early  times  in  the  detonation 

process  is  intended  to  represent  the  first  in  a  series  of  efforts  to  gain 

knowledge  of  the  basic  physical  characteristics  of  the  detonation-shock 

wave  formation  process.  The  primary  objective  of  the  experiments  was  to 

examine  the  validity  of  Sachd’  scaling  procedures  under  the  assumption  that 

the  airshock  is  Independent  of  the  explosive  details.  The  density  of  the 

TNT  charges  was  varied  from  a  density  of  l .01  gm/cc  (loose  powder)  to 

hWS  gm/cc  (cast  TNT)  with  the  result  that,  at  least,  the  airshock  pressures 

were  Independent  of  charge  density  and  dependent  only  on  the  energy  > eleased 

24 

by  the  detonation.  Kingery  h_  prepared  a  report  that  Is  primarily  in¬ 
tended  to  cover  the  administrative  and  gross  techni  al  information  for  the 
forthcoming  "Operation  Snowball”,  However,  this  report  does  include  tabular 
and  graphical  data  of  predicted  overpressure,  impulse,  positive  duration, 
arrival  time,  and  dynamic  pressure  versus  distance  for  50C  ton  TNT  surface 
detonation,  as  extracted  from  the  data  of  references  id  and  !9. 
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Ths  report  hy  Flshsr25  prat  ant*  overprtsiure  and  Impulse 

measurement*  fro*  detonations  of  bar*  apharleal  east  TXT  charges  6.5  Inchaa 

In  diameter  weighing  approximately  7.65  lb.  Thaaa  measurement*  are  compared, 

with  axeallant  agreement,  to  the  experimental  peak  pressure-distance 

26 

measurements  by  Itonar  and  Blaakney  and  with  the  theoretical  pressure* 

27 

dletanc*  curve  for  caat  TNT  generated  by  Kirkwood  and  Brinkley  .  Since 
the  data  of  Stoner  and  Blaakney  wai  obtained  from  denotation  of  Fertolite 
chargee,  the  coaler laon  sade  by  Fisher  In  thle  report  involves  the  computation 
of  the  magnitude  of  the  factor  required  to  convert  from  a  TXT  to  Fentollte 
bate.  For  acalod  distance*  ranging  from  6.5  to  11.  Fisher  finds  that  the 
weight  of  a  TXT  charge  required  to  prahice  the  same  overpressure  at  the 
same  scaled  distance  a*  generated  by  a  Fentollte  charge  varies  from  1.0* 
to  1.25.  As  noted  previously.  Baker2  recosneods  a  magnitude  of  1.1  for 
this  conversion  factor.  In  addition,  Fisher  compares  the  peak  overpressure 
computed  from  knowledge  of  the  shock  velocity  using  the  Bankine-Hugonlot 
equation  for  Ideal  gases*  with  the  measured  peak  overpressures,  again  with 
reasonably  good  agreement. 


*  See  Fart  B  of  the  Appendix  arc t lor  this  report  for  s  summary 

discussion  of  the  kenklns-ttugontot  relationships. 
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Reference  2$  discusses  in  aavaral  sections  Cha  methods  employed  to 


determine  various  blast  wava  pa  raw a  fro*  a  100  Ion  IKt  detonation  in  August 
1961  at  Cha  Sufflald  experimental  Station  blast  tasting  site.  Unfortunately,  no 
raaultant  data  ara  givan  in  thia  report.  However,  tha  r a port  by  Kingary,  at  al** 
praaanta  overpressure,  Impulse,  duration,  arrival  turn,  and  dynamic  praaaura  data 
naaaurad  during  tha  saaw  taat.  Thaaa  data  ara  eonparod  to  pradlctad  curvaa 
with  vary  good  agreement,  tinea  no  ground  reflection  factors  ara  givan  in  thia 
raport,  tha  raportad  data  can  not  be  uaad  to  provide  a  eat  of  "fraa-oir"  blaat 
parameters.  However,  onca  aueh  a  aat  of  free*alr  data  ia  determined  from  other 
aourcea,  tha  data  raportad  by  Kingary  at  al  can  then  be  »«ad  to  determine  tha 
ground  raflaetion  coafficiant(a)  applicable,  at  laaat,  ever  tha  range  of  acolad 
diatancaa  <1.67  £  Z  &A6.4  ft/lb  and  typo  of  Urrain  covered  in  thia  100 
ton  detonation. 

)0 

The  data  to  ba  presented  by  Ruetanlk  and  Lewis  in  a  forthcoming 
publication  was  collected  during  a  aeries  of  tests  designed  to  evaluate  the 
effaces  of  a  blast  environment  on  moving  airfoils.  The  results  of  these  Casts 
wars  obtained  tram  a  eat  of  carefully  calibrated,  fast  response  measuring  and 
data  recording  Instrumentation.  A  scaled  distance  range  of  1*1)  - 

covered  by  this  series  of  sight  data  collection  runs.  Overpressure,  d.ratlon, 
and  arrival  time  were  obtained  from  these  tests,  in  addition  to  data  concerning 
tha  response  of  the  test  aerlas  of  airfoils. 
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The  report  by  Fisher  le  presented  here  prlaerlly  to  tneorporece  the 
reflection  coefficient  dete  presented  in  reference  31  in  this  etudy.  Fisher's 
results  lndlcete  thet  the  measured,  average  reflection  coefficients  varied  from 
1.S3  for  e  herd  cley  surface,  to  1.99  far  a  perfect  reflector  slsuleted  by  the 
'nterseetlon  of  shock  fronts  from  two  Identlcc)  spherical  eha*g*t  flieo  slswl- 
eeneously.  The  tebuleted  results  lr  reference  31  ere  presented  es  functions 
of  the  eceled  height  of  the  ehnrge  ebove  the  groi  -d.  The  overpressure,  impulf*. 
end  duretlon  date  given  In  reference  32  ere  presented  for  e  seeled  dlatenee 
rente  between  epproxlsetely  1.3  ft/lb1^  to  22  ft/lb1^.  Unfortuneuly,  the 
ambient  pressure  end  tempers ture  ere  not  given  In  this  reference.  As  noted 
previously,  without  these  a«M>lent  conditions  the  dete  can  not  bo  applied  to  the 
evaluation  program  discussed  in  this  report. 


‘f;  i 


The  following  tsble  auauriiti  ti<«  Information  an.)  hUet  wave 
parameter  data  to  be  found  in  the  data  reference  eourcee  that  have  been 
diecueeed  in  thla  section.  A  definition  of  the  symbols  used  In  this  table 
follovs  the  .able. 
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Definition 


Aablmt  Prtkiurt 


Amblant  IifmtEi 


Sgsiatl  »"**■ 


°k  (or  °A> 


Z  -  tL/V 


Scaled  Dlataact  ft/ 
DUtasca  h«  Chart*  Cantar  ft 
tffaatlva  hrifUt  Weight  *f  Ch*-f*  lb 
Nak  SIKa-On  Ovarpraatvre  pal 


N*k  lifUcttd  Onrsnum 


Peak  IlMa  Pool  live  Iapvlee 


Mi  Sofia*  tad  Paaitiva  Iapvlae 


pat-— ec 


PartUla  Valaclty 


Ap<t> 


Shock  Vll«tt]r 

Paaitiva  frratU*  *f  Ovarpcaeaac* 
Shock  Arrival  Ttaa 
Tina  Rlatary  *f  Ovarpraeaar* 
flat  Rlatory  af  PartUla  Valaaity 


TIm  Rlatary  iaoaly 


\  (or  Ib/ft  ) 


17 


B.  Evaluation  of  Aval labia  Data  Sou ice a 
1 .  Ganara  1  Consents 

In  gcnaral  any  reference  that  provides  Information  on  the 
blaat  parameters  associated  with  the  detonation  of  an  explosive  charge,  aa 
exemplified  by  the  data  aourcaa  dlacuaaad  in  Fart  A  of  thla  aactlon,  doaa  not 
provide  complete  data  unl«*«  the  ambient  praaaura  and  temperature  (or  equivalent 
esfclent  round  apacd)  cordltiona  indar  which  the  data  waa  collected  are  aptclfled. 
Without  knowledge  of  there  ambient  condltione  the  acalad  blaat  parameter 
ralatlonahlpa.  summarised  in  Part  A  of  the  Appendix  aactlon  of  thla  report,  can 
not  bo  computed.  Direct  conparlaon  of  the  bleat  data  from  varloua  aourcaa  can 
not  be  aecompliahad  without  flrat  reducing  all  data  to  a  conalatent  act  of 
environmental  condltione,  aa  provided  by  the  Sacha*  acallng  ralatlonahlpa.  Hence, 
of  all  the  data  aourcaa  dlacuaaad  In  Part  A  of  thla  aactlon,  only  thoae  that 
define  the  magnitude  of  theae  ambient  conditlona  will  prev**  •!»«•»••’  in  rh#  later 
eectlona  of  thla  report. 

In  addition,  many  of  the  data  aourcaa  conaulted  during 
preparation  of  thla  report  do  not  report  "free-eir"  blaat  parameter  data;  i.e., 
theae  referencea  provide  date  that  Incorporatea  an  unknown  (since  It  la  unreported) 
ground  reflection  coefficient.  Such  data  ere  not  adequate  for  uae  In  defining  a 
consistent  sot  of  fre«-slr  blest  parameter  data;  however,  they  ere  useful  In 
establishing  the  magnitudes  of  the  ground  reflection  coefficients  once  the 
free*alr  blast  characteristics  have  boon  defined  from  other  sources  that  provide 
direct,  free-elr  date.  Of  those  data  sources  that  p.wvlde  free-alr  blast 
characteristics.  It  has  bean  noted  that  the  majority  of  the  experimental  over¬ 
pressure  end  shock  arrival  time  data  extracted  from  the  various  source  referencea 
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tend  to  agree  with  one  another  quite  well.  Unfortunately,  this  la  not  the  case 
with  reapect  to  poaitlve  phase  duration  time  data.  (The  date  of  Ruetenlk  and 
Lewia^0  la  probably  the  beat  aource  of  Internally  conslatent  duration  time  data, 
at  leaat  over  the  acaled  dlatance  range  from  Z  ■  5  ft/ lb  ^  to  Z  ■  15  ft/ lb  ^). 
Poaitlve  phaae  duration  time  date  la  the  critical  information  n*«dcd  to  generate 
a  aat  of  aalf-conalatent  bleat  cheracterlatlca,  alnce  the  veaultent  character  of 
the  particle  velocity  and  density  tlne-hlatorlea  will  be  dependent,  to  e  large 
extent,  upon  the  correlation  between  the  theoretically  computed  end  experimentally 
meaeured  overpressure  time  history  characteristics. 

Finally,  since  a  sal (•conslatent  set  of  blest  parameters 
resulting  from  Pentollte  detonations  he  a  been  presented  by  baker  and  Schuaman^, 
the  prlsmry  emphasis  of  this  report  hea  been  directed  toward  providing  e  similar 
set  of  parasasters  for  TNT  detonations.  Hence,  the  remainder  of  this  sectlom  will 
be  llmltad  to  e  specific  eveluatlon  of  the  TNT  data  aourcea  summarised  la  Table 
T,  Part  A  of  this  section. 


I 
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2.  Evaluation  of  Specific  Data  Sourcea  for  TNT  Deton'tiona 

The  experimental  data  reported  by  Welbull  ha a  been  uaed  In 
the  letter  sections  of  this  report  to  provide  e  aegmant  of  the  experlmentel ly 
determined  blast  fie!'*  parameters.  It  is  to  be  noted  that  the  Renklne-Hugonlot 
condition*  and. the  dete  of  Sheer  end  Day55  were  applied  to  Welbull *s  reported 
dete  in  order  to  determine  the  magnitude*  of  the  shock  front  overpressure  functions. 

Aa  previously  noted,  the  overpressure  dete  reported  by  Doll 
end  Salmon5  dies  not  agree  with  other  overpressure  date.  This  la  probably  due 
to  the  limiting  response-time  characteristics  of  the  measuring  instrumentation 
employed  in  reference  S.  The  data  presented  in  reference  6  is  also  applied  in 
this  study,  since  it  is  in  reasonable  agreement  with  other  similar  dete.  It 
is  to  be  noted  that  an  ambient  pressure  of  14.7  psl  was  assumed  for  purposes  of 
seeling  these  dete. 

Although  the  dete  sheets  presented  in  reference  8  do  not 
illustrate  actual  experlmentel  dete  points,  tha  curves  drawn  in  this  reference 
ware  used  in  the  preliminary  stages  of  plotting  the  final  curve*  given  in  this 
report.  This  was  done  merely  to  insure  that  the  general  trend  of  the  experlawntal 
data,  ami  of  the  theoretical  calculations,  did  not  deviate  by  large  orders  of 
Mgnitude  from  the  predictions  given  in  reference  8. 

The  various  curves  presented  in  reference  1)  ere,  as  noted 
previously,  based  on  the  experlmentel  date  reported  by  Welbull*  end  Fisher  end 
Fittman*.  Hence,  Implicit  use  of  the  date  contained  in  reference  1)  has  been 

*  The  Ranklne-Hugoniot  conditions  ere  valid  only  if  <  strength  of  the  shock 
is  not  sufficient  to  Induce  dissociation.  Above  overpressures  of  approximate¬ 
ly  4.5  e  taro  spheres ,  the  data  of  reference  31  which  accounts  for  dissociation 
effect*  should  be  used. 
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accomplished  In  t'-is  study.  In  addition,  the  duration  data  of  Grenstrom 

Illustrated  In  reference  13  has  alee  been  applied  In  thla  study. 

The  date  reported  by  Groves**,  and  amplified  It.  reference  17 
has  also  been  found  co  be  applicable  In  this  study.  However,  the  seme  duration 
date  reported  by  Dewey**  appear  too  high  by  approximately  a  factor  of  tuo.  in 
addition,  the  ground  reflection  coefficients  In  reference  16  do  not  agree  with 
those  of  reference  IS  (or  equivalently,  reference  17).  The  date  reported  In 
reference  IS  is  not  applicable  to  tha  study  discussad  in  this  report.  Tha  lata 
given  in  reference  19  is  for  a  ground  burst  detonation.  Unfortututaly,  no  ground 
reflection  coefficients  ere  defined  in  referance  19;  hanca,  these  data  can  not 
be  used  to  define  tha  frea*elr  characteristics  of  a  TNI  detonation.  However, 
once  the  experimental  ovarpressure-scaled  distance  ralatlonshlp  had  been  established 
from  the  data  of  references  6,  17,  and  23,  reflection  coefficients  were  computed 
to  effect  agreement  between  the  plotted  date  end  that  of  refarance  19.  Tha 


duration  end  arrival  tine  data  reference  19  were  then  modified  by  these 
reflection  factors  end  applied  to  the  appropriate  “free-eir"  time-seeled  distance 
plots  in  order  to  sx>re  completely  define  these  blast  parameters. 

As  noted  previously,  the  limited  dots  given  by  Dewey***  is  not 
in  a  fora  that  lends  Itself  to  appllcet'on  to  the  manner  of  presenting  the  deta 
given  in  this  report.  The  date  Illustrated  by  Sheer  end  Wright**  has  been 
obtained  from  references  6  end  19,  end  has  therefore  foramd  e  part  of  the  deta 
uaed  in  this  study. 

The  data  presented  by  Rudlir.**’**  Is  outside  the  region  of 
Interest  in  this  study  (reference  22)  end  can  not  be  satisfactorily  extracted 
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irom  reference  23  because  ui  inCeipulaCicm  difficulties  that  arise  duu  to  the 

24 

manner  in  which  the  data  are  reported.  The  turves  illustrated  by  Kingery 

are  not  applicable  to  determination  of  the  frae-ait  blast  anvlronrae.it,  since  no 

ground  reflection  coefficient*  are  defined  in  thia  reference,  the  overprassur* 

25 

data  reported  by  Kishet  is  applicable  to  this  study  uni  has  been  used.  The 
data  given  in  referenca  29  is  tor  s  ground  burst  and  doas  not  presant  ground 
reflection  coefficient  magnitude*.  These  date  were  used,  however,  in  the  same 

30 

manner  us  previously  discussed  for  reference  19.  The  data  of  Ruetenlk  end  Lewie 
are  accompanied  by  reflection  coefficient  magnitude*  and  era  therefor*  directly 

* 

applicable  to  establishing  tha  frcr-alr  blast  charactaristic*  of  a  TNT  detonation  . 
The  overpressure  duration  time  data  establlalied  in  reference  30  do  not  exhibit 
aa  much  scatter  as  do  other  duration  data  obtained  from  other  aourcee.  Unfortun¬ 
ately,  the  data  of  reference  30  only  span  a  scaled  distance  range  from  approxi¬ 
mately  5  to  15  ft/lh'^;  hence,  thr  theoretical-experimental  comparisons  evalu¬ 
ated  in  this  report  must  rely  on  less  reliable  duration  data  outside  tha  scaled 
distance  range  covered  in  reference  30. 


The  two  Pfotul  I  tv  denotations  included  to  th»-.-  .hit.)  agree  quite  well  with 
the  curves  presented  in  Reference  I. 
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lU.  AIR  BUST  PROPERTIES  INFERRED  FROM  THEORY 

We  ui.ll  discuss  In  this  section  the  techniques  employed  by 
various  Investigators  for  analytical  prediction  of  blast  waves  from  con* 
veutional  explosives.  Because  little  work  of  this  nature  was  published 
during  the  decade  1947-1957,  we  will  first  briefly  consider  the  theoretical 
work  accomplished  prior  to  this  time.  We  will  next  consider  more  recent 
(and  mare  complete)  predictions  In  some  detail. 

A.  Theory  Prior  to  194/ 

Before  World  War  II,  essentially  no  theoretical  predictions 
hsd  been  made  of  the  characteristics  of  blsst  waves  generated  by  any  type 
of  explosive  source.  Because  of  the  need  for  such  predictions  during 
the  war,  intensive  research  was  carried  out  in  blast  theory  for  both  eon* 
ventlonal  and  nuclear  weapons,  In  more  or  less  parallel  efforts.  The 
results  of  the  early  calculations  for  blast  from  nuclear  weapons  (the  so-called 
Los  Alamos  "M  Problem")  are  summarised  by  Fuchs  In  reference  34.  Since  these 
result*  <i.i  n.,t  strictly  apply  to  the  subiect  of  this  report,  we  merely  mention 
them  here  ior  sake  of  completeness.  The  tsaiurlty  of  the  early  theoretical 
predictions  of  air  blast  frow  convent ional  explosives  were  su.de  by  Kirkwood 
and  Brinkley,  in  research  conducted  under  sponsorship  of  the  tational 
Defense  Research  Council  (NDRC).  Their  work  Is  reported  In  a  number  of 
NDRC  monthly  progress  reports,  and  In  several  susssary  reports*^' 

Hakino^  has  summarised  the  Kirkwood- Hr inkley  theory,  and  coaspared  It  with 
such  experimental  data  as  existed  at  the  lime  of  his  report  (1951). 


2* 


From  Hit-  Ktrkwooc1- Brink  ley  theory,  one  can  obtain  predictions  ci  peax 

overpressure ,  positive  impulse,  and  positive  chase  duration  for  tree*  llr 

hursts  at  sea  level  of  spherical  Pentollte  and  TNT,  from  the  surface  of  the 

1/3 

explosive  chtape  out  to  scaled  distances  of  .lightly  lest  than  30  ft/lb  . 
!>’  predictions  can  be  made  of  time  histories  of  various  blast  parameters, 
at  either  fixed  or  moving  points  in  space. 


/ 
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1.  Brodtk  Theoretic* 1  Analyses 

Brode^'  has  developed  *  numerical  technique  for  the 

solution  of  the  flow  p*  remoter*  behind  shocks  vruch  Is  applicable  to  * 

number  of  different  esses.  In  getters  1,  the  technique  involves  the  ns*  of 

Ranklne-Hugonlot  condition*  in  e  numerics!  Integration  of  finite  difference 

equations.  The  difference  equations  sr*  solved  step  by  step  end,  in  ttie 

limit  of  lrfini teslmal  differences,  converge  to  the  exact  differential 

equation*  of  hydrodynamic  flow  in  Lagrsnglan  form.  To  evold  the  difficulties 

which  srlse  from  the  discontinuities  in  the  flow  parameters  when  traversing 

s  shock,  Brode  employes  the  artificial  viscosity  technique  developed  by  von 
42 

Neumann  srd  Klchtmeyer  .  The  computational  method  esn  bs  applied  ta  sny 
problem  Involving  the  spherical  flow  of  gasses  behind  s  shock  wave  If  one 
can  formulate  the  equation  of  state  for  the  gasses  Involved,  and  can  determine 
the  total  blast  energy  released  on  initiation  of  the  shock. 

A  detailed  discussion  of  the  problems  Brode  has  solved 
utiiiH  *.!•«  above  IvJinlque  is  presvnted  below. 

a.  Point  Source  Explosion  in  *•*  Ideal  Cae  ()  »  1_,4>^  ■ 

A  point  source  explosion  is  on*  produced  by  the 
Instantaneous  release  of  energy  from  sn  Infinitely  small,  concentrated 
source.  It  differs  from  the  blsst  generated  by  th  detonation  of  sn  explosive 
charge  In  that  no  high  temperature,  high  pressure  gases  are  generated,  snd 
th*  flow  consists  only  of  the  propagation  of  a  spherical  shuck  wav*  into 


25 


V 


► 


tin-  surrounding  medium.  Brodc  has  calculated  this  prcblet.  cut  to  a  shock 
overpressure  of  O.i  atmosphere  for  an  ideal  gas  with  /  ■*  1,4.  The  results 

presented  are  peak  overpressure,  peak  dynamic  pressure  and  peak  particle 
velocity  as  functions  of  shock  radius;  total  pressure,  particle  velocity, 
and  density  as  functions  of  lagranfar.  distance;  overpressure,  particle 
velocity,  density  and  compression  in  units  of  their  peak  values  as  functions 
of  tulerian  distance;  duration  of  positive  phases  for  pressure  and  particle 
velocity  and  positive  and  negative  impulse  at  functions  of  Eulerlan  distance, 
and  peak  dynamic  pressure  as  a  function  ot  time  in  units  of  positive  duration 

of  velocity.  The  significant  feature  of  these  results  la  their  departure 

43 

from  the  Taylor  strong  shock  curvet  at  overpressures  below  20  atmospheres. 

39 

b.  Point  Source  Explosion  in  a  Heal  &;s 

The  solution  to  this  problem  is  the  same  as  that 
for  the  ideal  gas  case  except  that  the  calculator  employs  the  equation  of 
state  for  a  real  gas  (air).  Computed  values  of  the  flow  parameters  are  pre¬ 
sented  in  the  asms  format  as  before.  Here  the  peak  overpressure-distance 
curve  lies  below  the  same  curve  for  the  ideal  gas  case  for  shock  pressures 

over  10  atmospheres.  Thereafter  the  iueal  gas  assumption  Is  reasonably  valid. 

40 

c.  Blast  from  an  Isothermal  Cas  Sphere 

The  gas  dynamics  arising  from  the  sudden  rele.se  of 
an  Initially  static,  high  pressure,  isothermal  gas  sphere  are  tundamentally 
different  from  those  for  a  point  source,  however,  the  general  computational 
method  described  before  is  also  applicable  in  this  cas.-.  Rrode  has  carried 
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out  Lit*"  coiupi>i.«t.lurt*  for  an  ideal  gas  at  three  different  initial  conditions 
ind  comp,  red  Llie  results  with  the  point  source  solution.  Two  spheres  were 
initially  at  2000  end  121  atmospheres  end  high  temperature  (norm*'  density) 
while  the  third,  at  121  atmospheres  Initial  pressure,  had  equal  Internal  and 
external  temperatures  (high  Initial  density).  In  f.Ms  case  tha  lssding 
•hock  it  strengthened  es  much  as  20  percent  after  a  period  of  time  dua  to 
the  effect  of  a  sacond  shock  which  overtakes  the  main  shock.  At  Increasing 
distances  from  the  origin,  tha  blast  wava  from  an  Isothermal  gat  aphera 
begins  to  converge  to  the  point  tourca  wave  and  the  two  will  colnclda 
(within  10  paresnt)  after  tha  shock  wava  haa  passed  through  a  mass  of  gat  10 
tlmss  at  largs  at  tha  Initial  mass  In  tha  tphara.  ?rtor  to  this  tha  ahock 
strength  Is  last  than  that  of  a  point  aourca  shock. 

d.  Blast  from  a  Spherical  Chaige  of  TNT* 

In  thla  case  the  numerical  calculation  ueas  tha 

44 

equation  of  etate  for  TNT  given  by  Jonaa  and  Miller,  and  an  aquation  of 

stetr  for  .ilr  determined  by  fitting  *  curve  to  the  computed  dais  of  Clllmnre 

and  HHsenreth  end  Becket.^  Initial  conditions  for  the  problem  are  the 

47 

detonation  front  conditions,  es  calculated  by  C.  I.  Taylor,  at  the  in¬ 
stant  the  front  retches  the  snrfece  of  e  bare  spherical  charge  of  TNT 
having  a  loading  density  of  1,5  The  air  outside  the  bleat  iroit  ie 

taken  to  be  at  standard  aee  level  condition*. 
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The  result i  nf  the  computations  yield  values  for  the  corvnon 
blast  va  e  parameters  both  at  the  shock  front  and  In  the  region  behind  the  shock, 
extending  back  into  the  negative  pressure  phase.  These  resuits  are  presented 
graphically,  and  for  peak  overpressure  versus  distance  they  are  comnured  with 
the  theoretic-i  results  from  Brodes  earlier  works  fo-  a  point  source  explosion 
In  air  and  In  an  Ideal  gas  of  ?  »  1.4,  and  with  Cranstrom's^  omplrical  curvt 
(or  TNT  at  1.S2  gm/cc  loading  danslty. 

It  la  specifically  noted  that  tha  calculated  resulta  art 
valid  only  for  bare  spherical  TNT  of  1,5  gr.t/ce  loading  density  detoneted  In  e 
standard  sea  level  atmosphere,  and  cannot  be  scaled  to  other  atmospheric  condi- 
tlons  without  a  simultaneous  and  proportional  change  In  the  charge  loading  denelty 
and  detonation  uavc  conditions.  In  addition,  the  results  seem  to  leek  internal 
consistency  as  evidenced  by  values  of  the  'ame  parameter  that  differ  by  more 
tluin  101  when  reed  (rom  different  curves.  S'nre  the  pe.per  does  not  repo~t 
results  In  tabular  form,  but  only  by  swans  of  graphs,  it  Is  not  possible  Co 
dctmrilnc  whether  this  Inconsistency  Is  inherent  in  the  ruin. latinos,  or  is  duo 
to  careless  plotting. 

e.  Spherical  Shock  Tub#  blasts 

The  numerical  calculations  are  carried  out  for  the  explosion. 
In  sir,  of  glass  filled  spheres  of  sir  and  helium  at  litgli  pressure  and  d.-nsity 
(ambient  temperature).  The  tumpucacions  and  results  are  analogous  to  those  of 
the  Isothermal  sphere  problem,  but  for  real  gasses,  and  are  mors  characteristic 
of  a  high  explosive  blast  than  an  Ideal  isothermal  gas  •.•here. 
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2.  Computations  of  M/kino  snd  Shear 
49 

Maklno  snd  Shear  have  developsd  a  computational  technique  to 
determine  free  air  blast  parameters  behind  the  shock  front  based  on  known  shock 
front  parameters.  The  method  involves  the  replacement  of  the  hydrodynamic  flow 
equations  for  the  spherical  flow  of  a  non- conductive ,  non-vlscous  Compressible 
fluid  by  s  system  of  ehsiacteristic  equations.  The  characteristic  aquations 
are  written  in  finite  difference  form,  using  first  differences  only,  and  are 
than  solved  numerically.  The  Initial  conditions  for  tha  computation  ara  the 
fraa  air  shock  front  parameters  for  spherical  50/50  Pa;>tollLs  darlvad  from 
Goodman's*  fitted  experimental  eurva  and  tha  Hugonlot  aquations.  1  For  scalad 
distances  greater  than  Z  ■  19. V,  tha  Klrkwood-Brinklay^  Asymptotic  solution 
Is  usad. 


The  authors  hava  limited  thalr  computations  to  tha  flow  region 
bounded  by  tha  shock  front  and  the  alr-axploslon  gaa  Intar taca,  due  to  uncartan- 
tlas  about  the  equation  of  state  for  the  axploalon  gasas.  For  acalcd  distance! 
greater  than  Zfe  2.04,  tha  computations  ara  only  carried  back  to  a  given  charac¬ 
teristic  curve  due  to  computer  limitations. 

Th '  results  are  presented  In  tabular  form  for  the  ranee  of 
acalad  distances  0. 132)9  Z  >  D.8  and  in  graphical  form  for  0. 1)2) IzS  1085. 
However,  the  authors  feel  that  the  calculated  reaulta  baaed  on  the  Kirkwood- 
Brinkley  asymptotic  expansion  are  questionable,  attributing  tha  error  to  Inherent 
deficiencies  in  tha  asymptotic  formula 

3,  Naval  Ordnance  Laboratory  Ccnaputationa 

The  Air-Ground  Explosions  Division  of  the  U.  S.  Naval  Ordnance 
Laboratory  has  recently  developed  s  computational  technique  to  calculate  all  the 
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blnst  wave  parameters  for  1  !h  of  TNT  detonated  in  free  air  at  sea  level 
SO 

conditions. 

'Jiu'  most  Important  now  approach  In  this  technique  la  the  nae 

2  J 

of  the  Umtau-Stanyukovitch  equation  ol  state  for  the  detonation  product*  of  TNT. 
Tills  equation  yield*  the  correct  Chapman* Jouguet  condition*  o:ul  also  tend*  to  an 
Ideal  a»s  with  1.24  a*  the  density  decreases,  lor  the  region  outside  the 
gone  containing  the  explosion  gaaea,  the  program  oaea  an  appropriate  equation  of 
state  for  real  air,  baaed  on  the  work  of  Hllsenrath.^ 

The  program  solve*  the  one  dimensional  ..ydrodyiutmic  equation* 
ot  spherically  symsetric  flow  using  standard  finite  difference  techniques.  (Options 
are  also  included  for  problems  involving  plane  and  cylindrical  symmetry.)  Shock* 
are  handled  by  the  artificial  viscosity  method. 

The  calculations  are  ear  rid  out  over  the  region  ext.n-^ng 
from  the  leading  stock  wove  back  to  the  charge  canter.  Thla  region  is  divided 
Into  300  xones  which  may  be  further  divided  among  a*  many  a*  30  different  region* 
or  materials.  A  variable,  internally  calculated  close  step  is  also  included. 

In  order  to  reduce  fomenting  1 11"  .  the  program  Includes  tour  type*  of  ramming. 

The  output  subroutine  print*  cable*  of  total  pressure, 
specific  volusns,  partlele  velocity,  temperature,  etc,,  versus  distance  at  fixed 
times.  Units  are  in  the  eg*  system  with  pressure  In  megehar*.  The  output  can 
also  he  presented  graphically  if  desired. 

The  subsequent  sections  of  this  study  will  use  this  compote* 

Clonal  technique  as  the  theoretical  basis  for  the  comparison  between  theoretlcsl 
predictions  and  experimental  meaaureoents  of  th.  Mast  •  >rameter*  from  free  air 


detona'.iuna  of  TNT. 
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IV.  COMPARISONS  OF  EXPERIMENT  AND  THEORY 


A.  Shock  Front  NriMtiri 

Tho  fro*  air  blest  ohock  front  p*-amot*rs  computed  by  the  MOL 
•Vundy  75"  Code  described  in  Section  III.B.3.  were  compared  with  compiled 
experimental  deco  extrected  from  the  vertOue  eourcei  dlseueeed  in  SecricM  il.A 
■no  II. 8.  A  detailed  diicueeion  of  the  reeulti  of  this  comparison,  end  the 
correction  factors  epplled  to  etak*  the  computed  velues  fit  the  experiment*! 
dots,  it  pretented  in  the  following  paragraphs. 

Sine*  the  artificial  viacoaity  method  used  in  the  theoretical 
calculation  of  the  thock  front  parameters  t»-ds  to  obteur*  tie  true  potltiom 
of  tha  shock  front,  it  was  assumed  that  the  instantaneous  time- specs  position 
of  the  shock  fron'.  coincided  with  eh*  peak  tabu  It  r  value  of  ovsrprsssurs  glvan 
in  the  computer  prine-out  sheets  (to*  Figure  1).  The  thock  arrival  time  deter¬ 
mined  by  this  approach  urns  then  plotted  against  sc»lsd  diseases  on  log- log 
paper,  Experimentally  determined  values  of  the  arrival  tioa  wets  next  plotted 
on  the  asms  sheet  in  order  to  determine  how  well  eh*  theoretically  predicted 
values  fit  the  experimental  data.  Tha  reauUi  show  that  tho  theory  la  in 
excellent  agreement  with  experiment  when  predicting  chock  errlvel  time,  and  no 
furehor  correction  ic  required. 

This  same  correlation  technique  war  chon  applied  to  tho  other 

•hock  front  parameter*  ,  ~  and  -t  ),  In  general.  Aha  theoretically 

r  *o  p* 

© 

predicted  magnitudes  of  these  parameters  were  significantly  lower  then  experimen¬ 
tal  values  at  tha  same  scaled  distances,  requiring  the  application  of  edditionel 
correction  factors  in  order  to  bring  the  theory  into  agreement  with  experiment. 
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T!>l>  necessary  correction  factors  for  each  parameter  were  determined  by  computing 
the  experimental  to  theoretical  ratio  fur  that  parameter  at  discrete  points 

ovu”  the  range  ot  scaled  distances  involved. 

,  AP  exp  I 

For  example:  C  A  ..  ■  s  I 

AV^or-  I  * 

In  general,  the  correction  factors  determined  by  this  approach  were  essentially 
constant  over  rather  large  ranges  of  scaled  distance.  The  veluee  determined 
for  the  correction  factors  to  be  applied  to  the  various  shock  front  parametars 
are  given  in  Tabic  11. 


Table  il 

blast  Wave  Parameter  Correction  Factors 
Parameter  C  Scaled  diatance  where  applicable 


Ap. 

i 

1.20 

Z  -  1 

Z  *•! 

0.916 

Z  se  l 

u 

a 

1.078 

ISZSl.l 

1.165 

5.1  -sZ 

1.089 

Z  25  5 

i 

a 

1  *#21 

,1.146 

It  should  be  noted  that  the  so-called  evpen mental  value*  of  peak 
particle  velocity  and  peak  density  were  calculated  from  the  experimental  data 
for  peak  overpressure  using  the  tables  given  hy  $hear  end  Uright^.  The  shock 
velocity,  which  tats  nor  calculated  by  the  NOI.  code,  was  determined  from  the 
corrected  peak  overpressure  using  the  same  nl>l. 


J2 
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B.  Time  Histories  ot  Hiaat  Parameters 

The  only  experimental  data  available  that  will  provide  a  companion 
with  theoretically  computed  time  hlatorlea  of  the  varloua  blaat  parameters  la 
the  tine  duration  of  poeltlva  ovarpreaaure.  Thla  comparison  waa  accomplished  In 
a  manner  analcgous  to  that  us«d  In  part  A  of  this  .taction.  The  results  a  ho  wed 
that.  If  the  theoretically  predicted  duration  of  poeltlva  phaaa, at  a  given 
acaled  dlatance,  waa  data  rained  aa  tha  tlaa  ftoat  the  atart  of  rlaa  in  praaaura 
ahead  of  the  a  hock  to  the  ties  tha  praaaura  again  returned  to  the  anblant 
condition)  the  theory  and  experiment  were  In  aatlafactory  agreement  for  acaled 
dlatancaa  laaa  than  11.3  ft/lb*^.  At  larger  acaled  dlataneaa.  thla  aw t hod 
tended  to  over  aetinate  the  duration  of  poaltlva  ovarpraaaura ,  requiring  the 
application  of  an  additional  corractlen  factor  to  the  start  of  rlaa  tins. 

The  corractad  tlaa  historian  of  ovarprassure  ware  then  obtained 
In  tha  following  awnnar:  It  was  assumed  that  tha  correction  factors  data  rained 
In  part  A  of  thla  section  also  applied  to  ovcrpresauras  behind  th«  shock  front, 
in  order  to  plot  lines  of  constant  ovarprassure  (corractad  for  praseur*  only) 
on  a  Z-t  plot,  the  desired  value  of  ovarprassure  was  divided  by  the  appropriate 
correction  factor,  yielding  the  corraspondlng  uncorracted  value  to  be  found  In 
the  computer  output  sheets.  This  Is  schematically  shown  In  Figure  1.  The 
Z  and  t  values  corresponding  to  glvan  values  of  ovarpraaaura  were  thei  used  to 
generate  the  Unas  of  pressure* corrected  constant  ovarpraaaura  on  a  Z-t  plot. 

In  addition,  Unas  representing  the  first  stuck  and  start  of  rise  of  praaaura 
were  placed  on  the  sane  plot.  Corractad  tlaw  histories  of  overpressure  were 
then  obtained  by  graphically  shifting  the  start  of  rise  curve  In  the  Z  dliectlon 
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Co  coincide  with  Che  first  shock  curve,  and  then  shifting  all  curvet  behind  the 

first  shod  th»  same  amount  In  Z  at  the  start  of  rise  curve. 

The  corrected  time  histories  of  particle  velocity  <u/#  )  and  density 

o 

ratio  (('/  )  were  then  obtained  in  ihe  same  manner, 

o 

C.  Ground  Reflection  Factors 

The  following  table  lists  the  reflection  fsctor(a)  either  given 
In,  or  computed  from  the  data  of  the  specified  reference.  These  ground  reflec¬ 
tion  factors  are  the  coefficients  used  to  compute  the  effective  free-alr  cheep 
weight;  t.e.,  the  weight  of  charge  required  In  the  absence  of  s  (ground)  reflect- 
Ing  eurfsce  to  generate  a  blast  wave  etrength  (overpressure)  equivalent  to  that 
produced  by  a  charge  In  the  presence  of  a  reflecting  eurtace. 

Table  Ill 

Gruuttd  Reflection  Factors 

Reference  Type  of  gurfe-e  Reflection  Factor 


1 

•  •  •• 

1.80 

17 

m  m  m  m 

1.60  - 

19 

.... 

1.89 

29 

.... 

1.39 

31 

hard  clay 

1.83 

•• 

perfect  reflector 

1.99 

•1 

water 

1.91 

it  ie  to  be  noted  that  only  reference  31  clearly  specifies  the 

magnitude  of  the  reflection  coefficient  as  a  function  of  the  type  of  ground 

30 

reflection  eurface.  The  data  of  Ruetenik  also  presents  reflection  factors; 
however .  the  magnitude*  of  these  factor*  indicate  that  a  modified  definition  of 
thl*  coefficient  ha*  been  ueed  In  this  reference. 
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0.  TNT-Pentolite  Convertlon  Faccort 

Table  iv  n*eii  TNT  to  Pentollte  convert  Ion  fectora  extracted  fros 


the  aourcet  cited  In 

the  flrat  column.  There  converalon 

fectore  are  ala^ly  the 

weight  of  TNT  neceaaary  to  obtain  the  tame  nverpreaaure. 

obtained  fr or  the  detonation  of  1  lb  of  Pentolttc. 

Table  IV 

iNfPentollte  Converalon  Pactora 

et  e  given  dlatanre,  ea 

Reference 

t  (ft/ib1/3)  Si 

mvtrelew  factor 

1 

m  mm 

1.1 

25 

4.5 

1.09. 

M 

5.0 

1.13 

H 

•.0 

1.1S 

H 

7.0 

1.22 

M 

*.0 

1.25 

11. 0 


1.25 


A  (teller  cosy eiieon  between  tW  experimental  data  for  TXT  need 
in  thle  report  end  Coodsan'a*  date  for  50/50  Pentellte  yielded  the  TXT  to  Pee* 
toltte  converalon  fectore  Hated  In  Table  V. 

Table  V 

TNT.pentollte  Converalon  Factora*Thla  Report 


i 

Converalon  1'ector 

.4*0 

1.20 

,M7 

1.04 

1.31 

1.10 
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Table  V  (Continued) 


Converelon  factor 


2.12 

2.96 

4.21 

j.bd 

11.3 

21.3 
47.0 


1.14 

1.15 
1 .18 
1.21 
i.oe 

1.07 

1.24 


baaed  on  theae  raaulta.  the  correction  fetor  of  1.1  given  In  reierence  1 
appeara  to  be  a  reaaonable  average  value. 
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B 


PRESENTATION  OF  SELF-CONSISTENT  DATA 


Thin  section  presents  large  acala  log-log  plots  of  the  self-consistent 
bleat  data  generated  by  thla  study  for  one  pound  of  TNT  detonated  in  free  air. 

Figure  1  preaenta  the  shock  front  values  of  the  various  blast  paraaMters 
at  functions  of  scaled  distance  trow  the  charge  enter. 

Figures  4,  i,  and  b  present  space-time  plots  of  overpressure,  particle 
velocity  and  density,  in  that  order. 

These  figures  are  to  be  found  at  the  end  of  this  rsport. 

La send  for  Flsures  3.  4.  5  and  6 

Symbols: 


■  side -on  overpressure,  atmospheres 

■  shock  velocity,  dimensionless 

■  particle  velocity,  dimensionless 

•  density,  dimensionless 
-  scaled  distance,  ft/lb1^ 

...  ■  seeled  time,  ms/lb1^ 


I  • 

w  ■ 

a  • 


a  * 


I  ■  scaled  positive  duration,  me/lb 

distance,  ft 
Wight  of  TNI,  lb 

sonic  velocity,  ft/sec 

a 

o 


Tffi 


1/5 
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P  ■  pressure,  pul 
P 

P  *  Q 
14.7 

Subscripts: 

•  «  shock  front 
o  ■  ambient  atmosphere 

th«  line  labeled  interface  on  Figure*  4,  5  and  o  designate*  the 
boundary  eeparating  air  froa  explosion  gas. 
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VI.  DISCUSSION 

In  till*  report,  w*  h*v*  presented  a  self- consistent  aat  of  blast 
wave  properties  that  characterize  the  parameter*  associated  with  detonations  of 
TNT  explosive.  This  self-consistent  sat  of  blast  wava  properties  has  bean 
generated  from  a  marriage  of  the  experimental  data  collected  and  recorded  by 
various  experimenters  with  th*  theoretic# l  computational  technique  developed  by 
personnel  of  the  Air-Ground  Explosions  Division  of  the  U.  S.  Naval  Ordnance 
Laboratory.  Whlta  Oak,  Silver  Spring.  Maryland.  The  datalls  of  thla  computational 
cechnlqua  are  described  In  refaranca  50.  Tha  experimentally  determined  data 
appllad  in  th*  thaoratlcal-axperimental  comparison  study  described  In  thla 
report  have  been  extracted  primarily  from  raferencaa  4,  6,  17,  19,  25,  29,  and 
50.  Other  sourcaa  of  a.'tperlmantal  data  end  other  theoretical  analyses  have 
been  ravlawsd  and  dlscussad  in  thla  report  primarily  to  provide  an  overall, 
comprehensive  examination  of  the  present  atete-ef-the-art  with  roapoct  to  bleat 
wava  characteristics. 

The  graphical  format  used  for  prasantatlon  of  th*  aelf-conatstont  set 
of  TNT  blast  wave  properties  wss  chosen  primarily  to  correspond  to  tha  foraat 
appllad  by  lakar  and  Schumen1  for  a  similar  presentation  of  Pantollte  properties. 
In  addition,  th*  acala  ot  those  graphs  has  beau  selected  to  provide  thrae-plac* 
reading  accuracy  fur  Che  entlra  scaled  dlstanca  range  covered  In  these  charts. 

During  the  affort  expended  i»  this  study  that  mbs  concerned  with 
evaluating  the  veri««s  sourcaa  of  axperlmental  date.  It  wss  clearly  avldent 
that  the  overpressures  and  errtval  time*  measured  by  different  experimenters 
tended  to  agree  quit*  well.  Howsvar,  tha  duration  time  measurements  were  quits 


scattered,  even  for  data  extracted  fro*  e  single  reference.  (Aa  noted  previously, 
tlit!  duretion  data  obtained  from  reference  30  we  a  the  baat  of  aoch  date  found 


during  this  study).  It  la  auggeated  that  future  experlawntal  blast  wave  detona* 
tlun  programs  stake  an  attempt  to  overcome  the  deflclanclaa  now  associated  with 
duration  time  oeasuresteiiCa.  Moreover,  It  la  suggested  that  in  future  detonation 
teats  the  measurement  of  ambient  pressure  and  temperature  be  eaphaaleed  ea  being 
aa  important  as  the  collection  of  the  bleat  data  Itself.  Much  of  the  data 
found  during  this  atudy  could  not  be  used  due  to  the  complete  lack  of  such 
Inf'  rsMtlon. 

The  curvee  preeented  In  this  rsport  are  fait  to  represent  the  beat 
fit  of  the  MOL  computational  results  to  the  existing  experimental  data.  Ae 
euch,  these  curves  can  be  uaed  to  ploa  futurt  expert  sun  tel  test  programs,  the 
data  from  which  can  then  be  uaed  ta  verify  the  results  Illustrated  la  thia 
report.  Aa  more  such  data  becomes  available,  le  may  be  evident  that  the  )corvee 
preeented  In  this  report  require  adjustment  end/or  revision.  To  Insure  the 
continuing  usefulnosa  of  thie  reoort,  It  la  suggested  that  the  curvaa  presented 
la  thia  raport  ha  updated  ••  required  by  additional,  future  data.  In  addition. 

It  la  suggested  that  consideration  be  given  ta  conducting  a  program  similar  to 
the  one  discussed  In  thia  report  for  Pentollta  explosive.  Although  reference  1 
presents  similar  data  for  PentoHt*,  these  data  are  baaed  on  a  marriage  of 
grode'a*  theoretical  calculations  for  TNT  with  Goodman**2  experimental  data  for 
PentoUt*.  Since  the  theoretical  analyst**0  applied  In  this  raport  oppoara  to 
agree  with  experimental  dare  better  than  the  theory  odvancod  by  Brode*,  end  alace 
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this  «ns lysis  could  be  recomputed  tor  Pontollto  explosive,  It  would  be  worthwhile 
2 

to  compere  Goodman's  dete  to  this  theoretlcel  approach.  Re-evaluation  of  the 
data  presented  in  reference  1  would,  It  Is  believed,  .  .  ova  eoate  of  the  Incon¬ 
sistencies  apparent  in  the  results  of  reference  1.  In  effect,  thle  re-evaluation 
would  provide  e  better  basis  of  comparison  of  Partcllte  end  TNT  detonations 
then  Is  available  at  the  present  time.' 
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APPENDIX 


A.  Blast  Parameters  Based  on  Bachs'  Scaling  Laws 

The  ratloa  presented  In  this  section  define  .he  v«leva,.c 
blast  wave  parameters  formulated  from  analysis  of  the  general  law*  of 
similitude  originally  proposed  by  Sachs^  and  caper lnentally  verified  by 
various  experimenters  (e.g.  reference*  2,  6,  7,  9,  10,  12).  Since  there 
«e  various  references  that  dlacuas  the  detail*  of  these  scaling  analyses,  they 
will  not  he  presented  In  thia  section.  (The  Interested  reader  la  re  fared 
In  particular  to  the  paper  by  Sperraxaa32  for  derivation*  of  the  scaled 
blast  wave  parameters).  Rather,  thia  brief  section  will  be  United  to 
Hating  the  form  of  these  blast  wav*  parameter*  aa  obtained  from  application 
of  the  general  similitude  relationship*  of  Sacha  that  account  for  variations 
in  ambient  pressure,  temperature,  and  aise  of  charge, 
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in  cuts  listing,  in  addition  to  the  symbols  defined  in  Table  I 


of  Section  1IA  of  this  report,  the  following  symbols  are  defined: 
P 


p  * 


Ratio  of  Ambient  pressure 
14.7  to  that  at  sea-level 


_  a 

a  I 

*  1139 


Ratio  of  ambient  sound  velocity 
to  that  at  sea-level 


It  is  to  be  noted  that  in  both  thla  Hating  and  in  Table  I 
of  Section  II  A  the  symbol  W  is  used  to  denote  the  "effective  free-air 
charge  weight. '*  by  this  Is  iseant  that  for  detonations  close  to  'the  ground, 
the  actual  weight  of  the  charge  i»  to  be  mul tipi led  by  the  ground  reflection 
factor  to  account  for  the  enhancesient  of  the  strength  of  the  blast  due  to 
the  presence  of  a  ground  (reflection)  plane. 


SI 


¥ 


a,  Runkine  *  Hugoniot  Relationship* 

the  basic  re latlunshl ps  among  the  properties  of  a  blast 

wave  arc  derived  trom  the  Rankine-Hugoniot  conditions  based  on  the  con* 

yervar.  m  of  .suss,  energy,  and  momentum  at  the  shock  front.  These  conditions, 

together  with  the  equation  of  state  for  air,  permit  the  relationships  between 

shock  velocity,  particle  velocity,  overpressure,  dynamic  pressure,  and  the 

density  of  the  air  behind  the  shock  front  to  be  derived.  These  relationships 

will  be  presented  in  this  sretion  without  derivation,  In  order  to  provide 

the  reader  with  expressions  relating  the  various  biast  wave  parameters 

that  are  applicable  for,  at  least,  shock  front  velocity  *  to  *  ambient 

sound  speed  ratios  leas  than  two  (overpressures  less  than,  approximately,  4.3 

atmospheres)*.  Fo  •  those  Interested  In  the  derivation  and  discussion  of 

these  relationships,  reference  to  the  texts  edited  by  Glasstonc*^  and 

34 

authored  hy  Couran*  and  Frirdrleks  should  he  made.  The  following  nomenclature 
lisllnx  dr i ine«  the  symbols  used  m  the  equations  presented  in  this  section. 
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ambient  sound  speed 
(ahead  of  shock) 

ambient  pressure 
(ahead  of  shock) 

peak  dynamic  pressure 

particle  velocity 
(behind  shock) 

peak  o*r. 
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relleeteil  ov.  .j.rt-isure 

shock  front  velocity 

•tensity  IwhinO  .shock 

ambient  density 
(alien. I  of  shock) 


tee  reference  J3f..r  the  technique."  t  »•»»!  *••  t  In  rases  exceeding  this 

limit. 
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The  relevant  equations  relating  thee*  parameters  art: 


The  following  curve  hat  been  extracted  fro*  ClaaetoM  55 
(ri|«re  ].N1  and  ia  presented  tn  this  report  In  order  to  Uluatrato  the 
character  of  these  equations. 
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to.  At  A  a  AtfUTY/UKfYATION  NO  TICK! 

Qualified  requestors  sty  obtain  copies  of  this  report  from  DDC 

it.  sveen BMSMTAar  eetes 

It  trONtOMIMS  MIUTARV  ACTIVITY 

DeF  tie*  Atomic  Support  Agency 

V^shlngton,  P.  C.  20301 

It.  4HTMCT 


Thl«  report  presents  in  a  convent...  ,  graphical  format  4  ss J f-sone latent 
»et  of  blest  vave  properties  for  •  T'.<1  explosive  Summarise  of  the  contents 
of  varir>v,  sources  of  pertinent  experiments!  bets  sre  given  in  this  report, 
together  with  evaluations  of  these  date.  *n  order  t;  justify  »he  use  of >niy 
a  selected  number  of  these  eourc.i  » ‘.•>*ueslon»  vf  various  theoretical 
methods  of  predicting  blast  vave  charac varieties  are  also  given  in  this 
report.  (U) 
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SELF- CONSISTENT  BUST  WAVE  PARAMETERS 


L  ORIOINATWO  ACTIVmft  btM  *ka  hm  m>4  iMthi 
af  iKt  HteMMita,  |IMH  Buinwil  af  D» 

Imh  activity  ar  atbar  aryaalaatlaa  IwyaMi  wllaf)  1 1  til  ay 
l Ha  iif  aw. 

la.  REPORT  IECWETT  CLAMflCATtOHi  Ehar  tb*  aver, 
all  aacarlty  alaaaUlcatlaa  af  Uw  rcyart.  latlaata  atwthir 
"Raatrlcta4  Dale**  la  laclntcA  MaAla*  la  la  be  la  aaa>v  A 
aaca  with  aff  npilaia  aaaarMy  ratalatlaaa. 

la.  GROUP:  Aalaaiaila  fcmyaRai  la  ayacNIai  la  DaO  > 
•aallaa  ICIA,  lb  Ml  Atwea  Patera  lafaatWal  Mawaal.  faa 
Iba  (ran#  ana* at.  Alia,  aha  ayyUcablr,  ahaa  Uw  ayliaiMt 
awblaya  Hava  baaa  vra4  far  Ovaaf  1  W  Or  try  4  aa  ahbar- 
lead. 


nwnwcnoaw 

wa  Hama*  by  aaaartty  .tfaaaUlcatlaa, 


(I)  “Qaallflc4  naaaalaia  aa  Jl 
rcyrrt  fraai  I'OC.” 

M  "Pantyi  Fanwaam  i  I 
rayart  by  Dt  C  la  aal  a«U. 

(D  “U  R  Oaaataaaat  ayaac  . 
Ibla  rayart  4  Cacti  y  Irrtt 
aaara  tbtll  aaaaat  ibraa. 


141  **U.  b.  -Ml.' aaaaclat  ->• 
*aya«l  directly  boa.  TT  * 
ahall  racial!  tbreafb 


m  iaa  raya.4  baa  at  •• 

Varvtaaa.  Bipatpaa  af  Caanna,  . 
aha  Ibla  (act  e*d  a*'ar  tbayrtea.  It  in 
U  yUPPVJMSPT/  **  >«T*»  l’  .  •  a 


I  REPORT  TITLE:  Ewer  iba  caaylala  rayaw  »«,  ta  an  •*“>!  rayaai  Ibrettyb 

aayllal  Inter*.  TlUaa  ia  all  aaaaa  tbevia  ba  aaalaaitfladi 

II  a  aaatla^il  ibla  aaaaat  ba  nlaala<  adibah  alaaalAaa.  _  ■ 

iter.  abaci  I  Mia  cUaatllaaltaa  la  all  aeyttatr  la  laaalnll  C*>  "W  ;licbllna(iMi.<; 

laretdtMaly  hlte~lr«  Iba  title.  I»h4  D"iC  a  the  aball  «  .• 

4.  DESCRIPTIVE  WTtl  n  ryyrcyrlhc.  carer  iba  af  _ _ 

tayan.  a.».  iMertra.  yreyraea.  aamaety,  aaaa;:.  at  *laal.  m  iaa  raya.c  baa  «av  • 

Ol«a  «Ha  iMlaaHa  4*tn  <abaa  a  ayuK1.  raye.iUe  ya— * :  lirvtcii.  Dciww  af  Caaaw,a, 

•*•***  call  ibla  (acl  and  war  iba  yrtaa,  I 

B.  AUTHuwtk  Ealav  Iba  aaaW  al  —>ai(b>  a»  data  aa  1 1  aiaat  raiiwt  •*  awia  I 

ar  la  Iba  ray  art.  Earn  law  aaaa.  Ibal  aaac,  ad  1 41  a  tatUal,  lary  nalit 

U  adlMery,  abaw raab  aa4 braacb a< earner.  Tba aaaa  a(  .  _ _ 

tbaytlaalyal  aalbar  la  aa  abaalala  alabraa  nhinaaab  II  POMOrvi  injTl 

(ba  dhawct.iur*  ■•••  • 

4.  REPORT  DATti  Erder  tba  beta  af  he  tryerl  aa  day,  |»*,.«r  raaaarab  ««4  4a 

anil,  yret,  ar  aahb.  yaaa  R aarr  tbaa aaa  dear  ryyear  _  .  .... 

aa  Iba  raped,  tree  dele  al  yabiurlUrv  I1AF.TPACTI  Babrr  eta 

la  lOTAL  MIMISR  Ob  P  AO  lb  Tba  la*cl  yeyc  aaaa  u  ^aa  a  Iaa  ijyaar  ,*aaabar 

abaaM  Mlaa  wwal  yaalaagaa  yyacalaaa.  la.  raw  tc  2£7n  addtSael  .yaw  h 
aaaber  af  ycyee  reharaiay  tadenreor*  bce.ieabad. 

lb.  RIMER  OP  ILVIXSJCS*  Cher  Aa  144  aa  bar  al  Ma  blabta  yaahabh  hr 

ralrraaaaa  rba4  la  tba  wyaH.  ba  aaalaaa Mai  Each  r  rat  ra  % 

•a  CONTRACT  OR  ORANT  RIMER:  Nayyrar  Wa.  aabv  art  ladlaaUaa  al  ha  wl. wry 
iba  ayyllaabla  aaabar  af  Iba  aarnraal  ar  Real  a*  car  ablab  hwauaa  la  ha  yarayayb.  ;  -** 
i  tba  iryau  area  arinaA  Than  la  aa  tbahallaa  aa  •  . 

IAb,bl4  PROJECT  NtMEh  Bala,  iba  ayyrayrUr  ha  i-ttaalal  ha«tb  h  -v 

ifrcy  brywawa  icailRaallaa.  ywb  aa  ywjy  awbar.  cry  vokm,  lu,^, 

ta  ORMIRATOR**  REPORT  WMBBRnh  Ealar  Iba  aMb  'aaaa  aabtaaha  caialafUt  t 

aial  ray  an  tawbar  by  attab  Iba  barraral  adU  ba  Mav.mat  aalaaaal  aa  hai  aa  aaaartty  «  - . 

aad  aaaarallat  by  Iba  arm  ah  bn  aauahy.  Tbta  nab  it  a»h  Hare,  aaab  aa  uilyaial  awta  1 
ba  ahaa*  H  iha  rayaw.  yra|aat  a*4a  aaaa.  arayraybla  a. 

aa  miiaa  rranar  aia.aah«b  If  tba - bw  baa*  ararta  bm  »IH  ba  MlaarW  by 
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